Numerical renormalization group study of two-channel three-impurity triangular clusters 
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We study triangular clusters of three spin-1 /2 Kondo or Anderson impurities that are coupled to two conduc- 
tion leads. In the case of Kondo impurities, the model takes the form of an antiferromagnetic Heisenberg ring 
with Kondo-like exchange coupling to continuum electrons. We show that this model exhibits many types of 
the behavior found in various simpler one and two-impurity models, thereby enabling the study of crossovers 
between a number of Fermi-liquid (FL) and non-Fermi-liquid (NFL) fixed points. In particular, we explore a 
direct crossover between the two-impurity Kondo-model NFL fixed point and the two-channel Kondo-model 
NFL fixed point. We show that the concept of the two-stage Kondo effect applies even in the case when the 
first-stage Kondo state is of NFL type. In the case of Anderson impurities, we consider the transport proper- 
ties of three coupled quantum dots. This class of models includes as limiting cases the familiar serial double 
quantum dot and triple quantum dot nanostructures. By extracting the quasiparticle scattering phase shifts, we 
compute the low-temperature conductance as a function of the inter-impurity tunneling-coupling. We point out 
that due to the existence of exponentially low temperature scales, there is a parameter range where the stable 
"zero-temperature" fixed point is essentially never reached (not even in numerical renormalization group calcu- 
lations). The "zero-temperature" conductance is then of no interest and it may only be meaningful to compute 
the conductance at finite temperature. This illustrates the perils of studying the conductance in the ground state 
and considering thermal fluctuations only as a small correction. 

PACS numbers: 75.30.Hx, 71.10.Hf, 72.10.Fk, 72.15.Qm 



I. INTRODUCTION 

Quantum impurity models describe localized single impuri- 
ties or impurity clusters in interaction with conduction bands 
of itinerant electrons. They appear in several different con- 
texts in condensed matter physics: as models for dilute mag- 
netic impurities in metals''^'-', as models of semiconductor 
quantum dots and other nanostructures embedded between 
conduction leads'*-^ and as effective models within the dynam- 
ical mean-field theory of bulk correlated electron systems''. 
Quantum impurity models are often studied also for their own 
sake due to the fascinating and rich behavior that they exhibit. 
Advances in the computational resources and improved im- 
plementations of the numerical renormalization group (NRG) 
technique^''*''''i'^'"'i2'i3.i4,i5,i6,i7,i8,i9,2o ^^^^ possible accurate 

and detailed studies of increasingly complex quantum impu- 
rity models featuring several impurities and several conduc- 
tion channels^''^^'^3-^4.25,26.27.28,29^ 

NRG permits to calculate finite-size excitation spectra, 
thermodynamic properties (such as impurity contributions to 
the magnetic susceptibility and entropy), various correlation 
functions (in particular zero-frequency correlators) and con- 
ductance through nanostructures both at zero temperature and 
at finite temperatures^". Such detailed knowledge about the 
behavior of the system under study at different temperature 
scales can be used to establish phase diagrams which delin- 
eate the parameter ranges with characteristic properties. In 
more technical terms, NRG allows to determine the possi- 
ble fixed points of the model, their stability with respect to 
various perturbation and crossovers between the fixed points. 
Depending on the nature of the excitations, the fixed points 
may be classified as either Fermi-liquid (FL) or non-Fermi- 



liquid (NFL) fixed points^SilLS. The excitation spectra of FL 
fixed points can be mapped one-to-one to the spectra of free 
non-interacting fermions (electrons); they have a characteris- 
tic appearance of equally-spaced lowest lying excited states. 
The excitation spectra of NFL fixed points cannot be related 
to non-interacting fermionic systems, they are typically more 
complex and the lowest lying excited states are not equally 
spaced (in some cases, however, the spectra may be described 
in terms of real Majorana fermions with twisted boundary 
conditions; the excitation spectra may then be given in terms 

0ffraCti0ns33-34,35.36,37.38)^ 

Non-Fermi-liquid properties of strongly correlated materi- 
als and proposed corresponding theoretical models have at- 
tracted the interest of the condensed matter community due 
to the very uncommon situation where the behavior of the 
system is radically different from what might be expected 
from the nature of its elementary constituents. The sim- 
plest models where NFL behavior emerges are quantum im- 
purity models such as the two-channel single spin- 1/2 impu- 
rity Kondo model^2i^2d£. In this model, the conduction band 
electrons attempt to screen the impurity moment as in the 
conventional Kondo effect, but as there are two conduction 
channels they tend to overscreen a single spin- 1/2 impurity. 
Since a strong-coupling state with overscreening is not sta- 
ble the system ends up, instead, in a non-trivial intermediate- 
coupling NFL state^". A mesoscopic system which exhibits 
the two-channel Kondo effect has recently been experimen- 
tally demonstrated^. 

In this work, we apply the numerical renormalization group 
techniques to study a complex impurity model which consists 
of three spin-1/2 Kondo-model-like impurities coupled by ex- 
change interaction so as to form a Heisenberg ring. Two of the 
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impurities are furthermore coupled to two conduction bands. 
This system is particularly interesting in that it has a rich 
phase diagram which includes fixed points known from sim- 
pler quantum impurity models. In contrast to simpler models, 
however, the three-impurity model allows to study crossovers 
between various non- trivial fixed points. Spurred on by a 
number of experimental achievement a^^''*^i"^'*''*^i^^i"^^'"^^i'*^ , the 
interest in the transport properties of three-impurity models 
recendy intcn^ificd^ "-^'-^^-^^-^"^-^^-^^-^^-^^-^^-^"-^'-^^-^-^-^'*-^^-^^ We 
thus also study the conductance through a system of three 
coupled quantum dots described by a related model featur- 
ing three Anderson-model-like impurity electron levels inter- 
connected by tunneling coupling and connected to conduction 
bands by hybridization. 

In Section|lI]we present the model and briefly discuss prop- 
erties of known limiting cases. In Section Hill we discuss the 
fixed points that are expected to characterize the various pa- 
rameter ranges. Numerical results for the Kondo-like model 
are presented in Section |IV] In this section we demonstrate 
the possibility of a direct cross-over between the two-impurity 
Kondo model fixed point to the two-channel Kondo model 
fixed point. At the same time, this result also establishes the 
validity of the two-stage Kondo screening concept in the case 
where the first stage of screening results in a non-Fermi-liquid 
fixed point. In Section |V] we discuss the transport through a 
triangular triple quantum dot systems connected to two con- 
duction leads. We emphasize that the notion of the "zero- 
temperature" conductance is of limited utility in systems with 
exponentially low energy scales, since experiments are per- 
formed at finite temperatures. 

II. MODEL 




Figure 1 : (Color online) Two-channel three-impurity triangular clus- 
ter. 

The three-impurity Kondo model with two conduction 
channels, represented schematically in Fig.[Tl is described by 
the Hamiltonian H = + i/jmp + H^, where 

ye{L,fl},fe,o-e{Ta} 

i<j6{1.2,3} 

i?c = Jk {Si -sl + Ss- sr) . 

The two {v = L and v = R) conduction bands are assumed 
to have linear dispersion, = Dk, where the dimensionless 
wave-number k ranges from —1 to 1 so that the bandwidth is 



2D. The impurities are described by the spin-1/2 operators 
Si with 1 = 1,2,3 and s^, is the conduction-band spin density 
at the position of impurity 1 (3) for v = L {R), 




where N is the number of conduction-band states, c ' , and 
c^ka are creation and annihilation operators for electron in 
band ly with spin a e {Tii}' and cTaa' is the vector of 
Pauh matrices (cTqq/ , ct^q/ , ct^q/)- Jk is the antifeiTomag- 
netic Kondo exchange constant. Finally, Jy parameterize the 
inter-impurity exchange interaction. We reduce the parameter 
space by considering only models with left-right mirror sym- 
metry (parity), i.e. J12 = ^23- While the parity-breaking is 
a relevant perturbation at some of the fixed points''^, its effect 
will not be studied in much detail in this work. We parameter- 
ize the exchange constants by 

Ji3 = Josin(/3V2), (3) 
J12 = J23 = Jo cos i(3n/2) . (4) 

Parameter (3 thus parameterizes the asymmetry (ratio) be- 
tween the exchange coupling between the impurities 1 and 
3 in the upper arm of the triangle and the exchange coupling 
between the side-coupled impurity 2 and the impurity 1 (or 3). 
The special cases are: 

9/3—1, which corresponds to the two-impurity Kondo 
model (plus one totally decoupled impurity); 

• /3 — 0, which corresponds to a linear chain of three 
impurities; 

9/3 = 1/2 which corresponds to a symmetric antifer- 
romagnetic Heisenberg ring, which is a magnetically 
frustrated system with two degenerate doublets in the 
ground state. 

Parameter Jq sets the overall scale of the inter-impurity ex- 
change coupling. 

In the two-impurity Kondo model^^ (i.e. (3=1 limit), there 
is a double Kondo screening regime for low inter-impurity 
exchange interaction J = Jq and an inter-impurity-singlet 

regime for high j67. 69,70,71,72,73,74,75,76,77 ^ rpjjggg pjj^ggg ^j.g 

separated by a quantum phase transition^^iS at J = J21K — 
cTjp where T^'' is the Kondo temperature for a system of a 
single impurity coupled to a single conduction channel with 
the same J^ (the proportionality constant c is of order 1; of- 
ten quoted value is ~ 2.2 obtained in the early NRG studies^-, 
however this particular value is not universal and true c de- 
pends on the details of the model and on the chosen defi- 
nition of the Kondo temperature). Exactly at the transition 
point, the system has a non-Fermi-liquid ground state and ex- 
hibits quantum criticality. This state is, however, unstable and 
for J — J2JK 7^ the systems flows to either of the two 
possible FL fixed points. We denote this NFL fixed point 
by 2IK. Recently, the two-impurity models were reexamined 
and it was shown that this fixed point is robust with respect 
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to parity, particle-hole symmetry breaking and various other 
asymmetries^. 

The system of three irnpurities in series (i.e. (3 — Q limit) 
has been studied in Ref. 80. This system has a non-Fermi- 
liquid ground state of the same type as the two-channel Kondo 

model (2CK)34-40,50,67.8 1.82,83.84.85,86.87.88,89.90.91 p^j. j^^ jntg^. 

impurity exchange interaction J, the local moment screening 
occurs in two stages: at the higher Kondo temperature T^^^ 
the local moments on impurities 1 and 3 are screened, while 
the local moment on impurity 2 is screened at an exponentially 
reduced second Kondo temperature 80.92 p^j. \i[g\i the 
three spins first lock into an antiferromagnetic spin-chain at 
T ^ J and the collective spin- 1/2 undergoes Kondo screen- 
ing at some lower temperature T2CK which depends non- 
monotonically on J**"'^^. The low-temperature two-channel 
Kondo (2CK) fixed point is stable with respect to particle-hole 
symmetry breaking^'*, but it is unstable with respect to parity 
breaking4<ia8a25. 

For /3 = 0.5 and sufficiently large Jq, the three impurities 
behave as a frustrated antiferromagnet with two degenerate 
ground state doublets at temperatures on the scale of Jo. The 
symmetry is broken by the coupling to the leads as there are 
only two conduction channels. It should be noted that in the 
more symmetric case of three conduction channels frustration 
induces a new type of non-Fermi-liquid behaviori2i2i. 



III. EXPECTED REGIMES 

The behavior of the system at various temperatures and 
inter-impurity coupling strengths is governed by the proximity 
to one of the following fixed points (see schematic represen- 
tations in Fig.|2]i: 

a) three independent local moments, LM, 

b) inter-impurity singlet (plus a decoupled spin-1/2 local 
moment), S, 

c) Kondo screening with 7r/2 phase shifts (plus a decou- 
pled spin-1/2 local moment), DK, 

d) two-impurity Kondo model non-Fermi-liquid fixed 
point (plus a decoupled spin-1/2 local moment), 2IK, 

e) frustrated antiferromagnetic Heisenberg ring, FR, 

f) antiferromagnetic spin chain with S = 1/2, AFM (two 
different fixed points exist, depending on J12 > J13 or 

J12 < J13), 

g) two-channel spin-1/2 Kondo model non-Fermi-liquid 
fixed point, 2CK. 

It should be noted that we have restrained ourselves to the 
fixed points which occur in the mirror (parity) and particle- 
hole symmetric case. In generic model, some of these fixed 
points are extended into lines or planes of fixed points. In 
particular, there appears a plane of Fermi-liquid fixed points 
(plus a decoupled spin-1/2 local moment) parameterized by 



two continuous quantities (phase shifts in even and odd scat- 
tering channel), of which the fixed points S and DK are special 
cases. There is furthermore a new plane of Fermi-liquid fixed 
points that is also parameterized by the two phase shifts but 
this time there is no decoupled local moment. The first plane 
(PI) occurs for (3 — 1 and the second one (P2) for /3 7^ 1; for 
/3 < 1, the crossover from PI to P2 occurs by a Kondo-like 
screening of the nearly decoupled local moment. 



a) LM e) 




S f) 




Figure 2: (Color online) Schematic representations of the possible 
fixed points. 

In Fig. [3] we present a schematic "phase diagram" of the 
system in the (J12, J13) [or, equivalently, {Jo,(3)] plane at 
very low (T 0) temperatures. Along the f3 — I line we 
find the S and DK Fermi-liquid fixed points, separated by the 
2IK quantum phase transition. Along the (3 — line as well as 
in the rest of the plane, the system always ends up in the 2CK 
fixed point. This is a special property of the parity-symmetric 
case (in the fully generic case, only PI and P2 fixed points are 
stable). For T > max( J12, J13), the system is in the local 
moment LM regime, while for T < min( J12, J13) the sys- 
tem is described by FR or one of the two AFM fixed points, 
depending on the values of J12 and J13. In the region of the 
parameter plane in close vicinity of the 2IK fixed point, i.e. 
for Jo ^ J21K (3^1, the system first approaches the 2IK 
fixed point, then crosses over into the stable 2CK fixed point. 

We emphasize that the 2CK fixed point dominates the low- 
temperature phase diagram, i.e. it is the stable fixed point for 
all /3 7^ 1. The values of parameters Jq and (3 affect only 
the way in which this fixed point is approached. A similar 
reasoning as in the /? = case*^" can also be applied to the 
general model. At some low enough energy scale, the im- 
purity cluster and the nearby conduction band electrons ef- 
fectively form a spin-1/2 object. This object is very local- 
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ized in the large- Jq limit (when the three impurities align an- 
tiferromagnetically into a spin-1/2 state) or extended in the 
small- Jo limit (in which case the composite spin-1/2 object is 
formed by the spin-1/2 of the impurity 2 and two collective 
spin-singlet states each consisting of an impurity and electrons 
which screen its spin). In both cases the coupling of this spin- 
1/2 object to the rest of the system (conduction band electrons 
at still lower energies in each of the lead) is antiferromagnetic; 
for large Jq this is explicit, since -terms in the Hamiltonian 
describe an AFM exchange interaction, while for small Jq 
this is expected by analogy with the two-stage Kondo effect 
in side-coupled geometry"* '^'''*'^. Thus the effective model 
takes the form of a single-impurity two-channel Kondo model. 
Since we consider reflection-symmetric models in this work, 
the exchange coupling to both channels are the same, thus the 
effective model with flow to the stable 2CK NFL fixed point. 
Of course, in the absence of reflection symmetry, the system 
would eventually end up in a Fermi-liquid fixed point (in the 
P2 plane). 
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Figure 3: (Color online) Schematic phase diagram of the triangular 
three impurity Kondo model at low temperatures. 

The presence of all the enumerated fixed points and the pro- 
posed features of the phase diagram are fully supported by 
the results of the numerical renormalization group calculation 
which are detailed in the following sections. 



IV. NUMERICAL RESULTS 

Calculations have been performed using the "NRG Ljubl- 
jana" package"^^. The impurity model, Eq. ([T), is particle-hole 
symmetric; in fact, it has a larger SU(2) isospin symmetry 
of which the particle-hole transformation symmetry is merely 
a subgroup2^i^2ii£. We performed all calculations taking ex- 
plicitly into account spin SU(2), isospin SU(2) and mirror Z2 
symmetry groups^' ' 2.20,25.97,98,99 _ -y^g j^^yg wssqA the discretiza- 
tion scheme described in Ref. 99 with the discretization pa- 
rameter A = 4. Averaging over four values of the twist pa- 
rameter z has been used'". The NRG truncation cutoff was 
set at the cutoff energy of E'cutoff = Qwjv, where tj^r is the 
characteristic energy scale at the iV-th NRG iteration, or at 



most 4000 states (which corresponds to approximately 32000 
states taking into account the degeneracies). To prevent sys- 
tematic errors, care is taken not to truncate within a cluster of 
nearby almost degenerate states. In all calculations presented 
in this article, we have used pJx — 0.2 which corresponds to 
the Kondo temperature T^^^ « 0.003D in the single-impurity 
Kondo model with the same parameter Jk. (Here we use Wil- 
son's definition of the Kondo temperature, i.e. for S* = 1/2 
Kondo model one has kBTKX^K) I {gp-sf = 0.701 or 
^bTkX^) I {sP'b)^ = 0.103. This definition is commonly, 
albeit not exclusively, used in NRG literature.) 

We first ascertain the presence of the expected fixed points 
by calculating the finite-size excitation spectra. These spectra 
can be represented in the form of the "renormalization flow di- 
agrams", some of which are shown in Fig.H] Flow diagrams 
show the NRG eigenvalue spectrum in units of the characteris- 
tic energy (or, equivalently, temperature) scale wat oc A~^/^ 
as a function of the NRG iteration number N (see Refs. p7lfT2l 
and ,20). We join the points by lines for easier visuahzation 
and interpretation. The colors (shades of gray) correspond to 
different sets of quantum numbers for the total isospin (/), 
total spin (5) and parity (P). The system is said to be near 
some fixed point when the eigenvalues do not change much 
between successive iterations (i.e. the lines are horizontal), 
while crossovers correspond to transitions between such re- 
gions. 

For A = 4, the single-particle eigenvalues which can be 
combined to give the Fermi liquid fixed point eigenspectra of 
a single Wilson chain with no impurities aro^ii^ 



77* = 0.8589029, 3. 99452, A^A^..., iVodd; 

f]* = 1.983281, 7.999996, A^/2,A^/^..., iVeven, 



(5) 



for odd and even iteration number N , respectively. In the first 
row of Fig.|4]we show the NRG eigenvalue flow for /? = 1 for 
three values of Jq. For J^/ D = 0.0154, we have Jq < J21K 
and we expect a flow to the Fermi-liquid double-Kondo (DK) 
fixed point with (5^™" = S"^^ = 7r/2 quasiparticle scatter- 
ing phase shifts. The 7r/2 phase shifts in both channels im- 
ply that the excitation spectrum of even-length Wilson chains 
correspond to that of odd-length non-interacting chain. In- 
deed, the lowest excitation energies are 77Q , 2rjQ , 3?/q , . . . to 
high accuracy. The quantum numbers are, however, differ- 
ent from those of the related two-impurity problem due to 
the presence of a decoupled spin-1/2 impurity (which also 
implies an additional two-fold degeneracy of all levels). For 
Jq/D — 0.0159, we have Jq > •^IiK' which corresponds to 
the flow to the Fermi-liquid inter-impurity singlet (S) phase 
with = 5'^^ = phase shifts. The excitation spec- 

trum of even-length Wilson chains corresponds to that of 
even-length non-interacting chain with lowest excitation en- 
ergy f^Q. Finally, the unstable fixed point for TV = 10, . . ., 20 
at Jo ^ J21K ~ 0.0158L» is the 2IK NFL fixed point with 
the energy spectrum (after suitable rescaling) described by the 
fractions 3/8, 1/2, 7/8, 1, . . 

The second row in Fig.|4]shows the flow diagrams for con- 
stant [3 — 0.92 for a range of values of Jq. In all cases the sys- 
tem ends up in the same stable fixed point, which is the 2CK 
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Figure 4: (Color online) NRG eigenvalue renormalization flow diagrams for even iteration numbers N. The states are indexed by the total 
isospin, total spin and parity quantum numbers. 



NFL fixed point with energy spectrum (after suitable rescal- 
ing) 1/8, 1/2,5/8, 1, 1+ 1/8, . . ., as predicted by the bound- 
ary conformal field theory approach to the 2CK problem^. 
Note the similarities in the flow diagrams in the first and sec- 
ond line, especially at high and intermediate temperatures (up 
to ^ 20). At lower temperatures, the coupling to the impu- 
rity 2 eventually drives the system to the 2CK fixed point for 
any value of /3 ^ 1 and Jq . 

Basic information about the magnetic correlations within 
the three impurity cluster may be obtained by considering 
the spin-spin correlations Cy = {Si ■ Sj) at zero tempera- 
ture. These provide insight in the competition between the 
inter-impurity interactions Jij and the impurity-lead Kondo 
exchange interaction Jk- In Fig.|5]we plot C13 and C12 both as 
a function of the interaction strength Jq and as a function of 
the asymmetry ratio f3. Note that if the impurities were decou- 
pled from the conduction channels (i.e. for Jk / Jo = 0), the 
spin coiTelation at zero temperature would depend only on the 
parameter (3, not on Jq. The behavior of the decoupled clus- 
ter approximates the properties of the impurity system in the 
large Jq limit, see the Jq/D = 1 plots in Fig.]?]?. At /? = 1/2, 
the decoupled cluster has C^y symmetry, thus C12 — C13. The 
two degenerate spin-doublets at this point are 

|a,) = l/x/2(|T,^,i>-U,^,T», 

|fo.> = i/V6(U,T,a> + |t,i,<T)-2|T,T,^))- 

For an equal mixture of these two states, we expect C12 = 
C13 = —1/2. For /3 7^ 1/2 the degeneracy is lifted. For 
/3 > 1/2, the decoupled impurities are in state |a), which 
corresponds to a local singlet state between impurities 1 and 
3 (ci3 = —3/4), while the impurity 2 is decoupled (ci2 = 0). 
For (3 < 1/2, the decoupled impurities are in state \b), which 
corresponds to a rigid antiferromagnetic spin chain (ci2 — 
-1/2, ci3 = 1/4). 

In the full model {Jk 7^ 0), the degeneracy between the 
doublets is lifted even at /3 = 1/2 by the coupling to the chan- 
nels. Curiously, in the limit Jq 3> Jk , the correlations C12 
and ci3 are not only different, but they even have opposite 



signs. This is explained by Fig.|5}): the point where C12 = C13 
is shifted from the value of /3 — 1/2. We also note that at 
large Jq, the transition between \ aa-) and \ bc) impurity ground 
states becomes increasingly sharp, thus a small change in /3 
leads to an abrupt change in the spin correlations. The spin- 
spin correlations curves at (3 = 1/2 play the role of separa- 
trix between two different limiting regimes (dashed curves in 
Fig. [5^. In the other limit of Jq ^ Jk, the spin correlations 
tend to zero due to the magnetic screening by conduction elec- 
trons. It may also be noted that the "side-coupled" impurity 2 
is aligned antifeiTomagnetically with respect to the "directly- 
coupled" impurities 1 and 3 for any values of parameters Jo 
and (3^1. Since antiferromagnetic exchange is a relevant 
perturbation, this implies that the local moment on impurity 2 
will always be screened (except for (3 — 1). 

At /3 = 1, we have studied C13 = (^i • S3) at the 2IK 
critical point, Jq = J21K' which is expected to be equal to 
— 1/4 due to a degeneracy between one singlet and one triplet 
state*"^ '™. Two general remarks concerning NRG calculations 
are in order at this point. The first concerns the value of the 
discretization parameter A. While even very high values of A 
typically lead to results which are qualitatively correct (in the 
single-impurity Anderson model, one can obtain decent mag- 
netic susceptibility curves even at surprisingly high A = 40), 
some details depend crucially on taking the A ^ 1 limit. The 
position of the quantum phase transition as a function of J and 
the value of {Si ■ S3) in the two-impurity Kondo model are 
one such example. At A = 4, we find J2Ik/-^ ~ 0.01585 and 
ci3 « -0.255, while at A = 2 we find J2*ik/^ ~ 0.01497 
and ci3 w —0.253. The convergence to the expected value of 
C13 = —1/4 is thus relatively slow. The second remark con- 
cerns the averaging over the twist parameter z (the "z-trick"). 
We find that the precise value of the exchange parameter J 
at the critical point depends slightly on the value of z at con- 
stant A, the more so as A is increased. To accurately study 
the detailed properties of the model in the vicinity of critical 
points, the averaging over z is thus better to be avoided and 
A should be kept small. Finally, it may also be noted that c 
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Figure 5: (Color online) Spin-spin correlation functions between the impurities in the upper arm of the triangle {Si ■ S3) and between the 
side-coupled impurity and one of the impurities in the upper arm (Si ■ S2). Note that Si ■ S2 ~ S3 ■ S2 due to reflection symmetry. The 
correlations shown are those at T = 0; in the major part of the parameter space, the final values are established on the scale ~ min( J13, J12) 
which is typically much higher than the Kondo temperature scales. 



in J21K = '^^K^ is here approximately equal to 5.0, not 2.2 
(still using Wilson's definition of the Kondo temperature), see 
Section II. 

We now focus on the vicinity of the 2IK fixed point, i.e. 
to the region f3 ^ 1, Jq ^ J2IK- ^'^^ /? = 1, the system 
crosses over from LM to 2IK fixed point at the two-impurity 
Kondo temperature y^^^^ . Since this fixed point is not sta- 
ble, the system then crosses over at some lower temperature 
rjf^^^ oc ( J - ^2Ik) V^^^^^^ to either S or DK fixed point, 
depending on whether Jq > J2jj^ or Jq < J2ik- This can 
be observed in Fig. |6] (black curve, f3 — 1), where we plot 
two thermodynamic quantities, the impurity contribution to 
the magnetic susceptibility ksTximpiT) / {gii%) ~ i^cS, and 
the impurity contribution to the entropy Si,„p (T) /fcs = In v^s 
as a function of the temperature. Note that the first quantity 
may be interpreted as the effective impurity cluster magnetic 
moment fi^s, while the second one can be related to the effec- 
tive number of degrees of freedom of the cluster, Vcff- At the 
Kondo temperature T^^^\ two local moments are screened, 
while the third decoupled spin remains free; correspondingly, 
/icff is reduced from 3/4 to 1/4. In the two-impurity Kondo 
effect, the NFL fixed point is associated with 1/2 In 2 residual 
entropy, thus the entropy goes from the LM value of 3 In 2 to 
In 2 + 1/2 In 2. The residual entropy of 1/2 In 2 is released at 



, when the system crosses over to the FL stable fixed 
point. Since impurity 2 is completely decoupled at (3 ~ 1, 
there is a residual spin- 1/2 local moment with In 2 entropy. 

Keeping Jo constant and decreasing (3, we weakly couple 
the impurity 2 to the rest of the system, see Fig.|6l One effect 
is the increased crossover temperature y^^^^) . More impor- 
tantly, the local moment on the impurity 2 is now screened 
at some lower temperature T^^. This can be observed both 
in //off curves, where the effective moment is reduced from 
1/4 to 0, and in the effective entropy which is reduced by 
1/2 In 2, which is a characteristic value for the two-channel 
Kondo effect. In our mirror-symmetric model, there is no fur- 
ther crossover and the 2CK fixed point is stable. Furthermore, 
by additional calculations we have checked that 2CK is the 
stable fixed point throughout the phase diagram. Fig. [3j for 
any Jq and /3 ^ 1, as expected. 

By analogy with the two-stage Kondo effect in the case of 
two impurities that are coupled to a single conduction chan- 
nel in a side-coupled configuration '^-'"'''"^ , the two-channel 
Kondo effect due to the side-coupled impurity is expected to 
occur on a temperature scale that depends exponentially on 
the effective exchange coupling of this impurity to the other 
two impurities, Joff, i.e. we expect a function dependence of 
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Figure 6: (Color online) Impurity cluster contribution to tiie mag- 
netic susceptibility and entropy for a range of values of /3 at fixed 
Jo. 



the form 



(2) 



K 



ij^ exp 



1 



Pcff Jeff 



(7) 



where T^^^ and T^^^ are first-stage and second-stage Kondo 
temperatures, while poff can be interpreted as the width of the 
band of the effective local Fermi-liquid quasiparticles result- 
ing from the first-stage Kondo effect, which is proportional to 
T^x ■ present situation, Joff oc cos(/37r/2). From the re- 

sults of calculations of the magnetic susceptibility along a line 
of values of /3 (from /? = 0.9 to /3 = 0.94) at a fixed value of 
the exchange coupHng J^jD = 0.0158, we extracted T"^^ 
using the prescription /ioff(T|-*-^^) = 0.07. Linear regression 
(see Fig.|7]) then gives 



^2CK 



-0.76-1.35 



1 



cos(/37r/2) 



(8) 



The good agreement confirms our anticipation that Kondo 
screening in multiple stages can occur whenever additional 
impurities are (indirectly) weakly coupled to the continuum 
of electrons via other (directly coupled) impurities. This oc- 
curs even in the case of more uncommon types of the Kondo 
effect. 

For Jq/J? = 0.0158, the system is already in the Fermi- 
liquid fixed point as the temperature begins to decrease to- 
wards the second-stage Kondo screening. We find, however. 
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Figure 7: The Kondo cross-over temperature as a function of the 
parameter /3 for fixed Jo/ D = 0.0158. 



that the second stage Kondo screening occurs even when af- 
ter the first stage of the Kondo screening the system cannot 
be described in terms of the Fermi-liquid quasiparticles, but is 
rather in a non-Fermi Uquid regime. In other words, the con- 
tinuum of NFL excitations may also serve as an impurity bath 
in a (second-stage) Kondo effect. By fine-tuning Jq around 
the critical value of J21K fixed (3 — 0.95, we have managed 



to reduce below \ see Fig. [8] This results in a 

direct cross-over from 2IK to 2CK fixed-point on the temper- 
ature scale of T(7, see the flow diagram in Fig. |9l It appears 
that the cross-over temperature Tc cannot be reduced to lower 
values by further tuning of the parameter Jq, which implies 
that the phase diagram in the (T, J) plane corresponds to that 
sketched in Fig. [TO] The 2IK regime does not extend down 
to zero temperature as would be the case (3 — 1 and there 
is a single stable 2CK fixed point at T = 0. As the temper- 
ature is reduced, we can either pass from 2IK to 2CK via an 
intermediate-temperature Fermi-liquid phase DK or S (path a 
in Fig. (To]), or directly from 2IK to 2CK (path b). As (3 is re- 
duced, the region in (T, J) plane that is governed by the 2IK 
fixed point, becomes smaller and eventually disappears. 
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Figure 8: (Color online) Impurity cluster contribution to the entropy 
for a range of values of Jo at fixed (3. is driven down until it is 
made equal to rl''^^. 



It has been shown that in the single-impurity two-channel 
Kondo model, the Kondo temperature is a non-monotonic 
function of the Kondo exchange interaction^. In the weak- 
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Figure 9: NRG flow diagram in the case of the cross-over from 2IK 
to 2CK regimes. 
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Figure 10: Schematic phase diagram in the (T, J) plane for constant 
parameter (3. 
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Figure 1 1 : The cross-over (Kondo) temperatures at which 

the system approaches the two-channel Kondo model fixed point as 
a function of Jq and 13. The numbers shown are the integer parts 
of the decadic logarithms, | logj^Q(r^'^^^/D)|. The omitted values 
for intermediate values of Jo correspond to the situation where it is 
difficult to determine T^2ck) since the cross-over occurs at relatively 
high temperature. 



a regular strong-couping Kondo fixed point in which the side- 
coupled impurity forms a Kondo state with electrons in either 
left or right channel, depending on which of the two effec- 
tive exchange constants is larger The intermediate 2IK fixed 
point, on the other hand, would be affected only little. 



coupling g — pqJk 'ti I regime we have 

while in the strong-coupling g ^ 1 regime 

with 7 — 30/4622.. By analogy, we expect that the temper- 
ature at which our three impurity system ultimately crosses 
over to the 2CK fixed point will be a non-monotonous func- 
tion of Jq at any value of (3, with maximum values occurring 
in the parameter range where Jq ~ T^^ with T^^ being the 
Kondo temperature of a single spin- 1/2 Kondo impurity with 
pJk ~ 0.2. This is indeed the case. In Fig.[TT|we show an 
overview diagram of the cross-over scale as a function of [3 
and Jq. The numbers displayed are the integer parts of the 
decadic logarithm of the cross-over temperature. For very 
small and very large Jq, T2CK is lower than the lowest tem- 
perature in the NRG iteration, ~ lO^^^D. It is worth empha- 
sizing that there is a relatively large section of the parameter 
space where the cross-over occurs at relatively high tempera- 
ture. This region is a continuation to the triangular impurity 
configuration of the intermediate regime found in the linear 
three-impurity model discussed in Ref. 

We may note in conclusion to this section that the par- 
ity breaking destabilizes the 2CK fixed point, while it is a 
marginal perturbation for the 2IK fixed point. A small parity 
breaking of the form Jk,l 7^ Jk,r or J12 ^ J23 will thus lead 
to another crossover from the now unstable 2CK fixed point 
to a stable PL fixed point (in plane P2) which corresponds to 



V. TRIANGULAR TRIPLE QUANTUM DOT 

A triangular triple quantum dot (TQD) model consists 
of three Anderson-like impurities interconnected by electron 
hopping and coupled to the conduction bands via hybridiza- 
tion (see also Refs. [13 and [sob . The geometry is still that 
displayed in Fig. [l] Various physical properties of systems 
with this kind of lattice connectivity have been previously 
studied^^'^^. In this work we focus on the Hamiltonian 

H = H]^ + iJinip + He 

where H^, remains the same as in Eq. ([T]l, while the impurity 
and coupling Hamiltonians become 

ka 

+ ^E4fc.rf3.+H.c. 

fc<T 

Here U is the on-site electron-electron repulsion, are hop- 
ping matrix elements, is the spin-cr electron creation op- 
erator on site i, and rii = dl^dia- is the occupancy oper- 
ator on site i. Finally, V is the hybridization matrix element 
which is assumed to be a constant independent of k. The hy- 
bridization strength is then F = wpV^ with the conduction 
band density of states p = 1/(213). There are no terms of the 
form eiUi which are typically included to describe the effect 



9 



of gate voltages. We instead assumed that the occupancy is 
near half filling, as implied by {U/2){ni — 1)^ terms. The 
particle-hole symmetry is broken by any finite hopping which 
makes the lattice connectivity non-bipartite and it should be 
noted that the deviation from half-filling may become signif- 
icant if hopping matrix elements are large, tij ^ U . We are, 
however, more interested in the limit of small ty ^ U . If 
the model parameters are chosen so that = ^tlj/U and 
pJk = ST/ttU, models ([U and ( fTTT i have namely very sim- 
ilar properties as long as tij ^ U. The essential difference, 
however, is that in model ([Til l the electrons are allowed to tun- 
nel from left to right conduction leads. It is known that inter- 
channel charge transfer destabilizes both the two-impurity and 
the two-channel Kondo-model fixed points^MS, therefore we 
expect that the ground state of the TQD system will be Fermi- 
liquid for any choice of parameters. We parameterize the tun- 
neling matrix elements using a new quantity a as 

ii3 = tosin(Q:7r/2), ^^^^ 
ti2 ^ t23 ^ to cos(a7r/2). 

Hoppings tij can be related to exchange constants Jy only in 
the limit tij <C U : then the relation between a and (3 is given 
by sin(/37r/2) cx sin^(a7r/2). 

The zero-temperature conductance through the TQD can be 
related to the quasiparticle phase shifts a s^^'^^-^^-"^^i'"'^-"^^ 

G = Gosin2(<5-™-5°dd), (13) 

where Go — 2e^/h is the conductance quantum. The phase 
shifts can be easily extracted from the renormalization flow 
diagrams when the system reaches the low-temperature sta- 
ble Fermi-Uquid fixed point27.2B,29.73,77,i04,i05,i06,i07.i08^ Phase 

shifts may be constrained in the presence of the particle-hole 
symmetry of certain kinds^. For a = 1, model ( fTTT l is 
particle-hole symmetric, however the quasiparticle scattering 
phase shifts are not fixed to any particular value; the stable 
fixed point therefore belongs to the plane PI, see Sec.|lll] For 
generic a, i.e. a ^ 1 and a ^ 0, the lattice on which the 
Hamiltonian (fTTT l is defined is not bipartite, which immedi- 
ately precludes any kind of particle-hole symmetry; the stable 
fixed point must thus belong to the plane P2. Finally, the case 
a = was discussed in Refs. 28 80 92: the phase shift in odd 
channel is constrained to 7r/2, while there is zero phase shift 
in the even channel (this fixed point also belongs to the plane 
P2). 

Properties of the double quantum dot (DQD) systems (i.e. 
a — 1 limit) have been studied in a number of works using 
various methoda^4i'"3,i09,ii",iiiiii2-ii3,ii4,ii5,ii6,ii7^ These cal- 
culations show that the conductance of DQD goes to zero in 
the limit of small t (as the dots become decoupled) as well 
as in the limit of large t (as the electrons occupy the bonding 
molecular orbital); the conductance peaks at the unitary con- 
ductance limit for t = t* such that J^s cT^^ where Jeg = 
it'^/U. The difference of phase shifts A = (5°™" - 6°'^^^ 
varies continuously and smoothly from to vr at i goes from 
to CX3 (it should be noted that scattering phase shifts are de- 
fined modulo tt). The fact that conductance becomes unitary 
at some some value of t is a simple consequence of A going 



through 7r/2; what is less trivial, but not unexpected, is that 

this occurs when Jeff ~ cT^iP ■ 

The conductance of TQD in series, on the other hand, was 
found to be unitary for all values of t, which is a direct 
consequence of the phase shift pinning at 5°*^^ = Tr/2 and 

srcvon _ n27.28.80.92 
^q.p. L^^-^-^— . 

At zero temperature the variation of the phase shifts as a 
function of t^ and a might be expected to be continuous and 
smooth throughout the entire parameter plane. Based on the 
considerations of the variation of phase shifts as a function 
of t in serial DQD and TQD problems, we can make qual- 
itative predictions about the dependence of the conductance 
on to and a. For to ^ t* we expect the conductance to drop 
monotonously from Go to as a goes from to 1 . For to ^ t* 
the conductance is Gq for a = and ~ Go for a = 1; it is 
reasonable to assume that the conductance drops to zero for 
some intermediate value of a. Finally, for to ^ t* we expect 
both a conductance zero and a conductance peak for some in- 
termediate values of a. 

The results of a NRG calculation of the "zero-temperature" 
conductance are shown in Fig. [12] We have taken into ac- 
count the spin SU(2) symmetry and parity Z2 symmetry. In 
these computations, a large value of the discretization param- 
eter has been used, A = 8; up to 150 NRG iterations were 
performed, which corresponds to an extremely low tempera- 
ture of r ^ 10^^^ D. Nevertheless, in some parameter ranges 
the stable Fermi-liquid fixed point has not been reached at this 
temperature. In the shaded box in Fig. [12] the excitation spec- 
tra at the last NRG iteration for some values of to were not 
those of Fermi liquid fixed points. Furthermore, to the right 
of the shaded box, in the immediate vicinity of a = 1, spec- 
tra obtained at the last NRG iteration correspond to Fermi- 
liquid fixed point in the PI plane (see Sec. HiH i. while it is 
known that the system should eventually cross over to a stable 
Fermi-liquid fixed point in the P2 plane. These findings can 
be explained by the presence of the exponentially low energy 
scales when the impurity 2 is nearly decoupled; see Eqs. (O 
and ([8]) for the equivalent behavior in the Kondo-like model. 
As a goes to 1, the cross-over temperature scale Ta becomes 
arbitrarily low, thus the "zero-temperature" conductance is ex- 
perimentally irrelevant. 

Transport experiments are decidedly performed at some fi- 
nite temperature Tcxp- As a is increased, at some value ajj 
the crossover temperature Ta will suddenly decrease (expo- 
nentially) below Tcxp- For all practical purposes, the quantum 
dot 2 will then be effectively decoupled from the rest of the 
system and its local moment will not be screened. It should 
be emphasized that a^i is less than 1 ; at finite temperature, the 
local moment effectively decouples from the rest of the sys- 
tem at a = au even though electrons are still able to hop on 
the quantum dot. 

Interestingly, this behavior also implies that the conduc- 
tance can be abruptly changed by driving the crossover tem- 
perature below the experimental temperature by relatively 
small variation of the gate voltages that control q'^. To our 
knowledge, this type of abrupt conductance change has not 
yet been experimentally observed. 

In Fig. [T3] we show the low-temperature conductance 
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Figure 12: "Zero-temperature" linear conductance through the upper 
edge of the triangular cluster and the quasiparticle scattering phase 
shifts as a function of the interdot-tunneling-asymmetry factors a for 
a range of tunneling parameters to. In the shaded box, the spectrum 
at the last NRG iteration does not correspond to that of a Fermi-liquid 
fixed point. The critical t* corresponding to the maximal conduction 
in the DQD model is approximately t* = 0.028(7. 



through the triangular TQD for a range of a as a function 
of tQ. Particularly interesting is how the deviation from the 
unitary conductance in the linear TQD (a — limit) develops 
when hopping between the first and the third impurity is al- 
lowed (a ^ 0). The differences are most significant for large 
values of to; this coincides with the regime of tij ^ U where 
the quantum dots are no longer restrained to half-filling for 
a 0. The change in scattering phase shifts can thus be re- 
lated to a change in the occupancy (Friedel sum rule). For 
larger values of a, the conductance may become zero at some 
value of to, see the a — 0.75 plot in Fig. [13] For a — 1 
we recover the known results for the conductance of the DQD 
system, but we must keep in mind that for a < 1, these results 
correspond to a finite temperature, Tcxp ^ Ta- 

VI. CONCLUSION 

The very rich phase diagram of the three-impurity two- 
channel system makes this a very useful toy model to study 
the possible behavior of generic two-channel quantum impu- 



rity models. We find all the fixed points familiar from sim- 
pler imp urity models and we find interesting crossovers, su ch 
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Figure 13: "Zero-temperature" linear conductance through the upper 
edge of the triangular cluster as a function of the inter-dot tunneling 
to for a range of tunneling-asymmetry factors a. 



as that between the two-impurity and the two-channel Kondo 
model non-Fermi-liquid fixed points. 

In a related triangular triple quantum dot problem we 
have demonstrated that the presence of energy scales which 
are extremely (exponentially) low implies that the "zero- 
temperature" conductance is experimentally irrelevant. In the- 
oretical studies of impurity clusters it is thus imperative to 
consider the thermal effects and to determine the temperature 
scale, below which the the transport properties of the system 
are determined by its ground state. 

It would be interesting to extend these studies to three- 
impurity three-channel models. Triple quantum dot systems 
with three conduction leads can be easily manufactured today. 
In order to determine the transport properties of such nanos- 
tructures in all parameter regimes (in particular in the low tem- 
perature regime where correlation effects and Kondo physics 
play a central role), an unbiased method such as NRG is re- 
quired. NRG becomes numerically highly demanding method 
in the case of multi-channel problems due to the high degen- 
eracy of quantum states that need to be considered. While 
two-channel calculations are now performed routinely, only 
very few three-channel calculations have been reported in the 
literature so far^^-^'* "^. Fortunately, the transport properties 
at low temperatures (if the system is near a Fermi-liquid fixed 
point) can be extracted from the energy levels of the NRG 
eigenvalue flow alone, which requires far less computational 
resources than calculations of thermodynamic and dynamic 
(spectral) quantities. 



Acknowledgments 

The authors acknowledge support of the SRA under Grant 
No. Pl-0044. 



11 



' J. Kondo, Prog. Theor. Phys. 32, 37 (1964). 

^ G. Griiner and A. Zawadowski, Rep. Prog. Phys. 37, 1497 (1974). 

' A. C. Hewson, The Kondo Problem to Heavy-Fermions (Cam- 
bridge University Press, Cambridge, 1993). 
L. I. Glazman and M. E. Raikh, JETP Lett. 47, 452 (1988). 

' L. Kouwenhoven and L. Glazman, Physics World Jan, 33 (2001). 

* A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Rev. 
Mod. Phys. 68, 13 (1996). 

' K. G. Wilson, Rev. Mod. Phys. 47, 773 (1975). 

* H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. 
Rev. Lett. 35, 1101 (1975). 

' M. Yoshida, M. A. Whitaker, and L. N. Oliveira, Phys. Rev. B 
41, 9403 (1990). 

W. C. Oliveira and L. N. Oliveira, Phys. Rev. B 49, 1 1986 (1994). 
" K. Chen and C. Jayaprakash, Phys. Rev. B 52, 14436 (1995). 

H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. 
Rev. B 21, 1003 (1980). 

" R. Bulla, A. C. Hewson, and T. Pruschke, J. Phys.: Condens. 

Matter 10, 8365 (1998). 
" W. Hofstetter, Phys. Rev. Lett. 85, 1508 (2000). 

F. B. Anders and A. Schiller, Phys. Rev. Lett. 95, 196801 (2005). 

R. Peters, T. Pruschke, and F. B. Anders, Phys. Rev. B 74, 245 1 14 

(2006) . 

" A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 076402 

(2007) . 

R. Zitko and J. Bonca, Phys. Rev. B 73, 035332 (2006). 
" A. I. Toth and G. Zarand, Dynamical correlations in the spin-half 

two-channel kondo model, arxiv:0801.4272 (2008). 
^° R. Bulla, T. Costi, and T. Pruschke, Rev. Mod. Phys. 80, 395 

(2008) . 

2' R. Zitko and J. Bonca, Phys. Rev. B 74, 045312 (2006). 
R. Zitko and J. Bonca, Phys. Rev. B 76, 241305(R) (2007). 

B. C. Paul and K. Ingersent, Frustration-induced non-fermi- 
liquid behavior in a three-impurity kondo model, cond- 
mat/9607190 (1996). 

M. Ferrero, L. de Leo, P. Lecheminant, and M. Fabrizio, J. Phys.: 

Condens. Matter 19, 433201 (2007). 

R. Zitko and J. Bonca, Phys. Rev. B 74, 224411 (2006). 

W.-z. Wang, Phys. Rev. B 76, 115114 (2007). 

A. Oguri, Y. Nisikawa, and A. C. Hewson, J. Phys. Soc. Japan 

74, 2554 (2005). 

A. Oguri and A. C. Hewson, J. Phys. Soc. Japan 74, 988 (2005). 
2'' Y. Nisikawa and A. Oguri, Phys. Rev. B 73, 125108 (2006). 
P Nozieres and A. Blandin, J. Physique 41, 193 (1980). 

C. M. Varma, Z. Nussinov, and W. van Saarloos, Phys. Rep. 361, 
267 (2002). 

I. Affleck, J. Phys. Soc. Japan 74, 59 (2005). 
J. Ye, Phys. Rev B 56, 1316 (1997). 

J. M. Maldacena and A. W. W. Ludwig, Nucl. Phys. B 506, 565 

(1997). 

" J. Ye, Phys. Rev Lett. 77, 3224 (1996). 
J. Ye, Phys. Rev Lett. 79, 1385 (1997). 

A. M. Sengupta and A. Georges, Phys. Rev. B 49, 10020(R) 
(1994). 

J. Ye, Phys. Rev B 56, R489 (1997). 

D. L. Cox and M. Jarrell, J. Phys.: Condens. Matter 8, 9825 
(1996). 

D. L. Cox and A. Zawadowski, Adv. Phys. 47, 599 (1998). 
'^^ R. M. Potok, I. G. Rau, H. Shtrikman, Y. Oreg, and D. Goldhaber- 

Gordon, Nature 446, 167 (2007). 
'^^ A. Vidan, R. M. Westervelt, M. Stopa, M. Hanson, and A. C. 



Gossard, Appl. Phys. Lett. 85, 3602 (2004). 
L. Gaudreau, S. A. Studenikin, A. S. Sachrajda, P. Zawadzki, 
A. Kam, J. Lapointe, M. Korkusinski, and P. Hawrylak, Phys. 
Rev. Lett. 97, 036807 (2006). 

F. R. Waugh, M. J. Berry, D. J. Mar, R. M. Westervelt, K. L. 
Campman, and A. C. Gossard, Phys. Rev Lett. 75, 705 (1995). 
M. Korkusinski, 1. P. Gimenez, P. Hawrylak, L. Gaudreau, S. A. 
Studenikin, and A. S. Sachrajda, Phys. Rev. B 75, 1 15301 (2007). 
D. Schroer, A. D. Greentree, L. Gaudreau, K. Eberl, L. C. L. Hol- 
lenberg, J. R Kotthaus, and S. Ludwig, Phys. Rev. B 76, 075306 
(2007). 

T. Jamneala, V. Madhavan, and M. F. Cronmiie, Phys. Rev. Lett. 
87, 256804 (2001). 

C. A. Stafford and S. Das Sarma, Phys. Rev. Lett. 72, 3590 
(1994). 

'^^ G. Chen, G. Klimeck, S. Datta, G. Chen, and W. A. Goddard, 
Phys. Rev. B 50, 8035 (1994). 

T. Kuzmenko, K. Kikoin, and Y. Avishai, Europhys. Lett. 64, 218 

(2003) . 

G. Zarand, A. Brataas, and D. Goldhaber-Gordon, Solid State 
Comm. 126, 463 (2003). 

T. Kuzmenko, K. Kikoin, and Y. Avishai, Phys. Rev. Lett. 96, 
046601 (2006). 

A. Oguri, Y. Nisikawa, and Y. Tanaka, Kondo screening of a 
high-spin nagaoka state in a triangular quantum dot, cond- 
mat/0606313 (2006). 

M. E. Torio, K. Hallberg, and C. R. Proetto, Europhys. Lett. 73, 
239 (2006). 

" Z.-t. Jiang, Q.-f. Sun, and Y. Wang, Phys. Rev. B 72, 045332 

(2005). 

M. L. Ladron de Guevara and P. A. Orellana, Phys. Rev. B 73, 

205303 (2006). 

" Z.-T. Jiang and Q. feng Sun, J. Phys.: Condens. Matter 19, 
156213 (2007). 

K. Kikoin, Low. temp. phys. 33, 153 (2007). 

A. M. Lobos and A. A. Aligia, Phys. Rev B 74, 165417 (2006). 
*° Y. B. Kudasov and V. M. Uzdin, Phys. Rev. Lett. 89, 276802 

(2002). 

V. V. Savkin, A. N. Rubtsov, M. 1. Katsnelson, and A. I. Lichten- 
stein, Phys. Rev. Lett. 94, 026402 (2005). 

B. Lazarovits, P. Simon, G. Zarand, and L. Szunyogh, Phys. Rev. 
Lett. 95, 077202 (2005). 

K. Ingersent, A. W. W. Ludwig, and I. Affleck, Phys. Rev. Lett. 

95, 257204 (2005). 

A. A. Aligia, Phys. Rev. Lett. 96, 096804 (2006). 

C. A. Busser, A. Moreo, and E. Dagotto, Phys. Rev. B 70, 035402 

(2004) . 

T. Tanamoto and Y. Nishi, Phys. Rev. B 76, 155319 (2007). 

1. Affleck, A. W. W. Ludwig, and B. A. Jones, Phys. Rev. B 52, 

9528 (1995). 

C. Jayaprakash, H. R. Krishna-murthy, and J. W. Wilkins, Phys. 
Rev Lett. 47, 737 (1981). 

B. A. Jones and C. M. Varma, Phys. Rev. Lett. 58, 843 (1987). 
™ B. A. Jones, C. M. Varma, and J. W. Wilkins, Phys. Rev Lett. 61, 

125 (1988). 

^' B. A. Jones, B. G. Kothar, and A. J. Milhs, Phys. Rev B 39, 3415 
(1989). 

B. A. Jones and C. M. Varma, Phys. Rev. B 40, 324 (1989). 
" I. Affleck and A. W. W. Ludwig, Phys. Rev Lett. 68, 1046 (1992). 
'''^ G. Zarand, C.-H. Chung, R Simon, and M. Vojta, Phys. Rev. Lett. 

97, 166802 (2006). 



12 



0. Sakai and Y. Shimizu, J. Phys. Soc. Japan 61, 2333 (1992). 
V. L. Campo and L. N. Oliveira, Phys. Rev. B 70, 153401 (2004). 
L. Zhu and C. M. Varma, Coherence in the two kondo impurity 
problem, cond-mat/0607426 (2006). 
M. Vojta, Rep. Prog. Phys. 66, 2069 (2003). 
M. Vojta, Phil. Mag. 86, 1807 (2006). 
R. Zitko and J. Bonca, Phys. Rev. Lett. 98, 047203 (2007). 

1. Affleck, A. W. W. Ludwig, H.-B. Pang, and D. L. Cox, Phys. 
Rev. B 45, 7918 (1992). 

H. B. Pang and D. L. Cox, Phys. Rev. B 44, 9454 (1991). 

H. Pang, Phys. Rev. Lett. 73, 2736 (1994). 

A. W. W. Ludwig and I. Affleck, Phys. Rev. Lett. 67, 3160 (1991). 
A. M. Sengupta and Y. B. Kim, Phys. Rev. B 54, 14918 (1996). 
^ J. von Delft, G. Zarand, and M. Fabrizio, Phys. Rev. Lett. 81, 196 
(1998). 

G. Zarand, T. Costi, A. Jerez, and N. Andrei, Phys. Rev. B 65, 
134416 (2002). 

I. Affleck and A. W. W. Ludwig, Nucl. Phys. B 360, 641 (1991). 
P. Coleman and A. J. Schofield, Phys. Rev. Lett. 75, 2184 (1995). 
P. Coleman, L. B. loffe, and A. M. Tsvelik, Phys. Rev. B 52, 661 1 
(1995). 

" N. Andrei and A. Jerez, Phys. Rev. Lett. 74, 4507 (1995). 

R. Zitko, J. Bonca, A. Ramsak, and T. Rejec, Phys. Rev. B 73, 
153307 (2006). 

C. Kolf and J. Kroha, Phys. Rev. B 75, 045129 (2007). 
''^ H. Kusunose, K. Miyake, Y. Shimizu, and O. Sakai, Phys. Rev. 
Lett. 76, 271 (1996). 

M. Pustilnik, L. Borda, L. 1. Glazman, and J. von Delft, Phys. 

Rev. B 69, 115316 (2004)^ 

R. Zitko, Nrg ljubljana, \http : //nrgl jubl jana . i j s . si/| 
(2006). 

H. R. Krishna-murthy, J. W. Wilkins, and K. G. Wilson, Phys. 
Rev. B 21, 1044 (1980). 

J. B. Silva, W. L. C. Lima, W. C. Oliveira, J. L. N. Mello, L. N. 



99 
100 
101 



103 
104 

105 
106 
107 



111 
112 



114 
115 
116 



Oliveira, and J. W. Wilkins, Phys. Rev. Lett. 76, 275 (1996). 
V. L. Campo and L. N. Oliveira, Phys. Rev. B 72, 104432 (2005). 
J. Gan, Phys. Rev. B 51, 8287 (1995). 

P S. Comaglia and D. R. Grempel, Phys. Rev. B 71, 075305 

(2005) . 

M. Vojta, R. Bulla, and W. Hofstetter, Phys. Rev. B 65, 
1 40405 (R) (2002). 

A. Georges and Y. Meir, Phys. Rev. Lett. 82, 3508 (1999). 

M. Pustilnik and L. I. Glazman, Phys. Rev. Lett. 87, 216601 

(2001) . 

W. Hofstetter and G. Zarand, Phys. Rev. B 69, 235301 (2004). 
L. N. Oliveira and J. W. Wilkins, Phys. Rev. B 24, 4863 (1981). 
L. Borda, G. Zarand, W. Hofstetter, B. I. Halperin, and J. von 
Delft, Phys. Rev. Lett. 90, 026602 (2003). 
P. Mehta, N. Andrei, P. Coleman, L. Borda, and G. Zarand, Phys. 
Rev. B 72, 014430 (2005). 

T. Pohjola, J. Konig, M. M. Salomaa, J. Schmid, H. Schoeller, 

and G. Schon, Europhys. Lett. 40, 189 (1997). 

T. Aono, M. Eto, and K. Kawamura, J. Phys. Soc. Japan 67, 1860 

(1998). 

W. Izumida and O. Sakai, Phys. Rev. B 62, 10260 (2000). 

D. Boese, W. Hofstetter, and H. Schoeller, Phys. Rev. B 66, 

125315 (2002). 

R. Lopez, R. Aguado, and G. Platero, Phys. Rev. Lett. 89, 136802 

(2002) . 

R. Aguado and D. C. Langreth, Phys. Rev. B 67, 245307 (2003). 
W. Izumida and O. Sakai, J. Phys. Soc. Japan 74, 103 (2005). 
C. Karrasch, T. Enss, and V. Meden, Phys. Rev. B 73, 235337 

(2006) . 

J. Mravlje, A. Ramsak, and T. Rejec, Phys. Rev. B 73, 241305(R) 
(2006). 

L. De Leo and M. Fabrizio, Phys. Rev. Lett. 94, 236401 (2005). 



